The Annals of Statistics 

2009, Vol. 37, No. 6A, 3529-3554 

DOI: 10.1214/08-AOS610 

© Institute of Mathematical Statistics, 2009 



ON NONPARAMETRIC AND SEMIPARAMETRIC TESTING 
FOR MULTIVARIATE LINEAR TIME SERIES 
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University of Tokyo and Tohoku University 

We formulate nonparametric and semiparametric hypothesis test- 
ing of multivariate stationary linear time series in a unified fashion 
and propose new test statistics based on estimators of the spectral 
density matrix. The limiting distributions of these test statistics un- 
der null hypotheses are always normal distributions, and they can be 
implemented easily for practical use. If null hypotheses are false, as 
the sample size goes to infinity, they diverge to infinity and conse- 
quently are consistent tests for any alternative. The approach can be 
applied to various null hypotheses such as the independence between 
the component series, the equality of the autocovariance functions 
or the autocorrelation functions of the component series, the sepa- 
rability of the covariance matrix function and the time reversibility. 
Furthermore, a null hypothesis with a nonlinear constraint like the 
conditional independence between the two series can be tested in the 
same way. 

1. Introduction. One of the main purposes of multivariate stationary 
time series analysis is to elucidate intrinsic relationships between different 
component series. Frequently, these relationships can be expressed in terms 
of specific constraints imposed on the spectral density matrix. For instance, 
the spectral density matrix of a separable time series is a product of the 
contemporaneous covariance matrix of the component series and the scalar 
spectral density function, which is common to them (see, e.g., Haslett and 
Raftery [14], Martin [23], Cressie [6], Guyon [13], Shitan and Brockwell [36] 
and Matsuda and Yajima [25]). If the underlying time series is Gaussian, 
the independence between the component series is equivalent, so that the 
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spectral density matrix is diagonal for all frequencies (see, e.g., Wahba [40]). 
The conditional independence of two components, given the others, is equiv- 
alent, so that the corresponding partial spectral coherence is identical to zero 
(see, e.g., Dahlhaus [7]). The time reversibility is characterized by the spec- 
tral density matrix being real- valued (see, e.g., Chan, Ho and Tong [3]). 
Furthermore, the equality of the autocovariance functions or the autocorre- 
lation functions of the component series implies that their spectral density 
functions are equal to each other with an appropriate scale-shift for the 
latter case (see, e.g., Coates and Diggle [5] and Diggle and Fisher [10]). 

In this paper, we formulate nonparametric and semiparametric hypothesis 
testing on multivariate linear time series in a unified fashion and propose new 
test statistics based on estimators of the spectral density matrix. First, we 
construct two estimators of the spectral density matrix, of which the first is 
always consistent and the latter is consistent only when the null hypothesis is 
true, and, next, we introduce a function to measure the discrepancy between 
these estimators. Then, the discrepancy is asymptotically standard normally 
distributed under the null hypothesis as the sample size goes to infinity, 
whereas, if the null hypothesis is false, it diverges to infinity. Related ideas 
are applied for discriminant analysis of time series by Kakizawa, Shumway 
and Taniguchi [18] and for nonparametric testing on univariate time series 
by Taniguchi and Kondo [38] and multivariate time series by Taniguchi, 
Puri and Kondo [39], respectively (see Taniguchi and Kakizawa [37] for a 
comprehensive exposition). Hong and White [17] considered specification 
testing on nonparametric and semiparametric regression models in a similar 
way. 

On the other hand, for parametric models, Paparoditis [29] and [30] and 
Delgado, Hidalgo and Velasco [8] considered test statistics based on another 
discrepancy measure between the hypothesized spectral density function (or 
matrix) and their estimators. 

The advantages of our test statistics follows. First, our test statisc- 

tics can be applied to both nonparametric and semiparametric hypotheses 
on multivariate linear time series in a unified way. They are robust against 
wrong decisions caused by misspecification, which parametric approaches 
often suffer. 

Second, under the null hypothesis, the limiting distributions of the test 
statistics are always normal distributions. Moreover, since they are scale 
invariant, that is, independent of the unit of measurement for the observa- 
tions, the expectations and variances of the limiting distributions can be 
expressed in a simple form, interestingly, some of them are known constants 
in the preceding examples. Hence, they can be implemented easily for practi- 
cal use. In contrast, the limiting distributions of some test statistics, which 
have already been proposed for specific hypotheses mentioned above, are 
unknown or have rather complicated forms, particularly when a nonlinear 
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constraint is imposed on the null hypothesis, which makes it difficult to show 
the appropriate validity of their significance tests (see, e.g., Dahlhaus [7]). 

Finally, if the statistics are normalized to be standard normal asymptoti- 
cally under the null hypothesis, then, under any alternative, they diverge to 
infinity and consequently are consistent tests. However, a trade-off for this 
advantage is that our test can only detect local alternatives of 0{n~( l+ ^l 4: ) 
with l/2</3<3/4 being slightly slower than ra -1 / 2 , which is shared with 
the test for regression models proposed by Hong and White [17]. 

The paper is organized as follows. We introduce the mathematical formu- 
lation of the null hypothesis and the test statistics in Section 2. In Section 
3, we derive the theorems on limiting properties of the test statistics. In 
Section 4, we apply the theorems to the examples mentioned above. Some 
simulation results are shown in Section 5. Section 6 is devoted to prove the 
theorems. The technical lemmas are listed in Section 7. 

The mathematical details on the lemmas, the theorems, the examples and 
the computational simulations are available from the authors upon request. 

2. Null hypotheses and test statistics. Let Z t = (Z±t, . . . , Z rt )' be an r-di- 
mensional zero-mean stationary Gaussian time series possessing a spectral 
density matrix /(A). /(A) is defined as periodic with the period 2tt for 

(-7T,7T]. 

Throughout this paper, A a b and A ab are generic symbols for the (a, 6)th 
element of matrices A and A~ l , respectively. A' is the transposed matrix of 
A. 

Let g(0,y) = (g a b(8,y)) be an r x r matrix-valued function, where y = 
(yap) is an r x r matrix and 8 = (8\, . . . ,8 V )' is a v -dimensional vector of 
parameters. We assume that, under the null hypothesis, /(A) satisfies the 
equation 

(1) f(X)=g(8,f(X)) 
for all A G (— vr,7r]. 

If gab(8,y) = y a b, for all a and b, it implies that there is no constraint 
imposed on /(A). Otherwise, (1) introduces some relationship between the 
components of /(A). The hypothesis is semiparametric for v > and is non- 
parametric for v = 0, respectively. 

Now, we define the test statistics. Given n observations Zi,...,Z n , in- 
troduce the discrete Fourier transform W a {\) = J2t=i ^at exp(ftA) for 
a = 1, . . . , r. The cross periodogram of Z a t and Z^ t is 

Iz,abW=W a (X)mW, a,6=l,...,r, 

and the periodogram matrix is -Tz(A) = (Iz,ab(X)). Iz(X) is defined as peri- 
odic with the period 2tt for A ^ (— tt,tt]. Denote the Fourier frequency ^p. as 
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Xj for j = 0, ±, 1, For notational simplicity, W a (Xj), Iz,ab(Xj) and Iz(Xj) 

are denoted by W a ,j, Iz,ab,j and Iz,j, respectively. Actually, they depend on 
the sample size n but we suppress it. 

Put ft = /(Aj), and define the unrestricted estimator of ft by the smoothed 
periodogram matrix f U{t = (fu,ab,t), where 

l m/2 

(2) fu,t = — w j J z,t+j, i = l,...,[n/2], 

j=-m/2 

and Wj (j = —m/2, . . . ,m/2) is a weight sequence with w* = Y^j l =- m /2 w j 
and [x] is the integer part of x. 

Under the null hypothesis, we define the estimator f^ t = (fR, a b,t) by 

fR,t = 9(0nJu,t), t=l,...,[n/2], 

where 9 n = (6i n , . . . , 9 vn )' is an estimator of 6. 

We can expect that, under the null hypothesis (1), /c/,t//j| is close to 
I r , the r x r identity matrix, whereas they are far from each other if (1) is 
false. Hence, we introduce a function that measures the discrepancy between 
/c/,t/i?t an d Ir- Let K(A) be a nonnegative function defined on all r x r 
complex matrices that are similar to a positive definite matrix, and let it be 

zero if and only if A = I r . Note that fu,tfs\ ls s i muar to f r^J 2 fu,tf r^J 2 ■ 
Then, we introduce the statistic 

[n/2] 

t=i 

where M t = (m a6;t ) = fu,tf^ v 

Now, we give a few candidates for K(A) (see [18]). The first one is 
Kullback-Leibler (KL) discrimination information defined by 

K I {A)=tr(A)-logdet(A)-r 

(see Kullback and Leibler [20] and Kullback [19]). Next, the J divergence, 
which is a symmetric version of the KL discrimination information, is 

Kj{A)=K I {A) + K I (A- 1 ) 

(see Kullback and Leibler [20] and Kullback [19]). Finally, using the Chernoff 
information (see Chernoff [4] and Renyi [33]), Parzen [31] proposed 

K a (A) = logdet(aA + (1 - a)I r ) - alogdet(A), < a < 1. 

We shall show in the subsequent section that, under the null hypothesis, 
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is asymptotically standard normally distributed with some constants rj and 
a. Consequently if fj n and a n are consistent estimators of 77 and a, respec- 
tively, and f\ n — rj = o p (m 1 / 2 /n 1 / 2 ), then the test statistic 




is also asymptotically standard normally distributed. Actually, rj and a 1 of 
some nonparametric examples mentioned in Section 1 are known constants, 
and they require no estimation procedure. On the other hand, T n diverges 
to infinity if (1) is false and, hence, is a consistent test for any alternative. 



3. Theorems. Before we proceed to show the main theorems, we intro- 
duce some assumptions and notation. First, we introduce the following as- 
sumptions: 

(Al) Zt is an r-dimensional zero-mean stationary Gaussian process; 
(A2) /(A) is a positive definite matrix for all A G (— 7r,7r]; 
(A3) /(A) is twice continuously differentiable for all A G (— 7T, 7r]; 
(A4) m = 0(vr), \ < /3 < |, and the weight sequence Wj,j = —m/2, . . . , 
m/2 is 

Wj = u(j/m), j = -m/2,..., m/2, 

where u{x) is a positive continuously differentiable even function on [—1/2,1 /2] ; 

(A5) g(0,y) is three times continuously differentiable for 6 and y; 

(A6) i^(^4) is four times continuously differentiable for A. 

Next, put O = (#10, • • • Ao)', gW = (s(0o,/(A)) and & = £(A t ) if the 
null hypothesis (1) is true. Finally, let k^J, be the second partial derivative 

d 2 K(A) / dA a b dA C( i evaluated at A = I r . The third partial derivative k^ cd e ^ 
is defined similarly. 

By the results in Appendix A. 13 of Liitkepohl [22], for Kj(A), 

( o\ J2) = fl, iia = d,b = c, 

{ ] ab > cd \0, otherwise. 

For Kj{A) and K Q (A), their second partial derivatives are given by 2n ab cd 

and (a — a 2 )^^ crf , respectively. 

From now on, to avoid abuse of notation and give the theoretical results 
in a concise form, we consider Kj{A) and assume (3), because, for other 
K(Ays, the corresponding results are given in the same way without any 
modification. 

Then, our main results are the following. 
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Theorem 1. Under assumptions (A1)-(A6) and (3), if the null hy- 
pothesis (1) is true and 9 n — 9q = O p (n _1//2 ) ; then the limiting distribution 
of {m/n) l / 2 {T n — (n/m)rj) is N(0,o~ 2 ), where 



E 



a, f3, 7,v=l 



VT Jo 
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a,f3,~f,v=la',f3'rf>,v'=l u 



(A) 

x {g aa '{X)gp'p{X)g^'{\)g V 'u{X) 
+ 9a-j'{X)g v ' p{X)ha'{X)gp'u{X)) dX 



and 



Vaf3jvW = - tr 



+ 



g (X)^—(X)g(X)i^(X) 



dy a p 
dg- 1 



dy 



1 

+ 2 



dg 



dy a fs 



(A) 



dy. 



-(A) 



71/ 



/3a 



+ _/7 (A) ^ (A)) 



and C u and D u are defined in Lemmas 9 and 10, respectively. 

Hence, if the null hypothesis (1) is true and % — r\ = o p {m 1 / 2 /n 1 / 2 ) and 
a 2 — a 2 = o p (l) , then the limiting distribution of T n is the standard normal 
distribution. 



Theorem 2. Under assumptions (A1)-(A6) and (3), if the null hy- 
pothesis (1) is false, and 9 n , fj n and a n converge in probability to some 
constants 6* , rf and a*, respectively, so that g(6*,f(X)) is a positive defi- 
nite matrix for all X in (— tt,tt] and a* is positive, then, for any sequence 
{C n },C n = o((nm) 1 / 2 ), 

lim Pr[f n > C n ] = 1. 

n— >oo 

Theorems 1 and 2 assure that our test statistics are consistent for any 
alternative. Furthermore, if we put 

°* = 2^Io K (f( X )9(0*,f(X))~ l ) dX, 

the proof of Theorem 2 shows that, if the null hypothesis is false, T n /^nm 
converges in probability to C* as n — > oo . 

Next consider asymptotic behavior of the test statistics under local alter- 
natives. To state our theorem, we define a class of local alternatives 

ff«»:/n(A) = /(A) + — LrAA), 
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where /(A) satisfies (1) with 6 = 6q and /*(A) is a spectral density matrix. 

Then, we introduce the following assumptions, which specify the behaviors 
of /*(A) and the parametric estimator 9 n : 

(A7) /*(A) is a positive definite matrix and twice continuously differen- 
tiable for A on (— 7r,7r]; 

(A8) Under the local alternative H an , there exists a nonstochastic se- 
quence {#*} such that 6 n — 6>* = O p {n~ 1 / 2 ) and linin^oc^mn) 1//4 (#^ — #o) = 

Po = (pio, ■ ■ ■ ,Pvo)'- 

Before we proceed to the theorem, we introduce some notation. Set g n (X) = 
g{9* n ,f n (X)) and g tn = g n (\ t ). r] n and a 2 n are defined by substituting g n (X) 
for g{X) in r\ and a 2 of Theorem 1, respectively. Then, we have the following 
result. 

Theorem 3. Under (A1)-(A8) and (3), if H an is true, fj n —rj n = o^m 1 / 2 / 
n 1 / 2 ) and a 2 —a 2 = o p (l), then the limiting distribution of T n is N(^/a, 1), 
where 




dg(\) 



ddi 



a, /3=1 



dg(X) 
dy a p 



2-i 



x-l/ 



dX. 



Theorem 3 implies that, although our test is consistent for any fixed alter- 
native, a cost for this advantage is that it can only detect local alternatives 
of 0((nm)- 1 / 4 ) = 0(?i^( 1 +/ 3 )/ 4 ), with 1/2 < /3 < 3/4 being slightly slower 
than n~ 1 / 2 . A similar feature is shared with the tests for regression models 
proposed by Hong and White [17]. 

Remark 1. We make several comments on the results. 

(i) We make an explicit comparison of ( A4) , the bandwidth of smoothed 
periodograms and the rate of convergence with those of related works. By 

1 /2 

replacing w* with the integral m J_ 1 ^ 2 n(x) dx and rewriting u(j /m) as 
u(nXj/(2irm)), we obtain 

^ I 

W * mj^{ 2 /2 u(x)dx \2irmJ 

On the other hand, some authors use K(Xj/h)/(nh) (see, e.g., Paparoditis 
[29] and [30]) or M n K(M n Xj)/n (see, e.g., Eichler [11]) with akernel function 
K(x) instead of Wj/w* . h and M n are called a smoothed bandwidth and an 
effective number of frequencies, respectively. Then, we have the relation 

m ~ nh ~ n/M n 
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as n — > oo, and our \Jmj nT n corresponds to T n of [30] and Qt of [11], respec- 
tively. In terms of m, (A4) is stronger than the assumptions like < [5 < 1 
(see [29] and [30]) or 1/2 < (3 < 1 (see [11]), which is the cost for considering 
nonparametric and semiparametric testing hypotheses in a comprehensive 
way but not a specific one. 

(ii) There are other test statistics alternative to T n . First, in some appli- 
cations, it may be better to leave out the frequencies where the determinant 
of /(A) is near zero to make a test statistic more stable. It suggests modifying 
T n to 

[n/2] 

where 4>t = 4>(^t) an d <^>(A) is a nonnegative weight function (see also [18]). 
The limiting behaviors of T n ^ are obtained in the same way as T n . 
Second, we introduce the quadratic function 

T Q , n = W^[{M t -I r f] 
z t 

and define the test statistic Tq^ by 




A similar idea is proposed by Kakizawa, Shumway and Taniguchi [18] for a 
discriminant analysis of multivariate time series. 

Finally, breaking up the frequency axis with nonoverlapping blocks and 
using only M (t _ 1)(m+1)+m/2+1 (t = l,...,L), we define T* by 

L 

T n = ^2 K ( M {t-l)(m+l)+m/2+l), 
t=l 

where L = [n/2]/(m + 1). 

Then, the test statistic is given by 

where B u = (J_y 2 u(x) 2 dx) 2 / J_i/ 2 u(x) 4 dx. For instance, Wahba [40] ap- 
plies T* to Example 2 of Section 4 by substituting Kj(A) for K(A). 

Similar to the proofs of Theorems 1 and 2, it is shown that the limiting dis- 
tributions of TQ tn and T* are iV(0, 1) under the null hypothesis (1), whereas, 
if the null hypothesis is false, Tq n j \Jnm and T*/^/nm converges in prob- 
ability to C* Q = i^/J r tr[(/(A) 5 (r,/(A))- 1 - I r ) 2 ]d\ and y/2D u /B u C*, 
respectively. 
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(iii) For practical use, we need a reasonable criterion to choose a specific 
test statistic among K{A) mentioned above or their alternatives. Theoret- 
ically, various concepts of asymptotic relative efficiency (ARE) of one test 
relative to the other one have been proposed. They differ from each other 
in intuitive appeal, or the availability of mathematical tools and efficiency 
comparison among tests cannot be done in a single ARE (see Serfling [35] 
for a comprehensive survey of AREs). Hence, the choice of a test statistic 
among candidates, from a practical point of view, is left to future studies. 

Here, we shall give a remark. According to Pitman's ARE approach (see 
Pitman [32] and Noether [27]), we consider the asymptotic power function 
when local alternatives converge to the null hypothesis as the sample size 

(2) 

goes to infinity. Because Kj (A) , K j (A) and K a (A) have the same n ab cd up to 
the constants, it follows from Theorem 3 that the expectations £/a of their 
limiting distributions under local alternatives are identical to each other, 
which means that they have the same asymptotic efficiency in Pitman's ARE 
sense. Similarly, Tn sn has the same asymptotic efficiency, whereas T* is less 
efficient than these test statistics, because its expectation of the limiting 
distribution is \J2D U / 'B u £/ 'a and, by the Schwarz inequality, 2D U /B U < 1. 



4. Examples. First, we show some generic formulas, which are helpful 
for deriving r\ and a 2 in Theorem 1. Since 

(4) ^w = _ r . (A) |w r . (A) , 



we can evaluate /x a/ 3 7! ,(A) without calculating 



9y a i3 



■ By applying (4), 



(A) 



(5) 



1 

-tr 
2 



OVa/3 Oy lv 



r . (A )fW r . (A) 

OVa/3 



oyju 



1 



+ ~/3 7(A)r * (A) 

(3a * 



X]Ma/37^(A) J sa y- 



8=1 



Then, for the calculation of rj, we have 

1 r 



(6) 



Y] ~ Ma/? 7 i/,2(A)5™,(A)5 7/ g(A)£L\ 

= ~ /^/37^3(A)^(A)<?7/3(A)dA 

a,/3,-y,u=l 71 J ° 
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Jo 





dX 

a/3 



and 



1 r 2 
(7) V - / /i a/ a 7l /,4(A)3a„(A)5 7/ 3(A)dA: 

. a ... 1 JO 



a,/3,7,i/=l ^ J ° 

Hence, only / lJia^u,i{X)gau{X)g 1 p{X) dX may requires a laborious evaluation. 
A similar technique reduces the evaluation of a 2 to a simpler one. 

Now, we only consider the three examples of those mentioned in Section 
1 and proceed to calculate rj and a 2 of them. For simplicity, we assume that 
u(x) = 1 on [-1/2, 1/2] and Wj = l. Then, C u = 1/2 and D u = 1/3. 

Example 1 (A separable model). The spectral density matrix of a sep- 
arable stationary process is expressed in the form 

/(A) = E/(A), 

where £ is an r x r positive definite matrix and /(A) is a scalar-valued 
nonnegative integrable function in (— 7r,7r] (see Matsuda and Yajima [25]). 
Set v = r 2 and 6 = vec(E), where vec transforms an r x r matrix into an 
r 2 -dimensional vector by stacking the columns of the matrix underneath each 
other. If we define g(6,y) = l(Ea=i Vaa/vaa^, then /(A) =g(9,f(X)) = 
F(S=l/aaWKa)S = S/(A). 

Next, we observe that g 1 (9,y) = r(YJ a=1 y aa / a aa ) *E 1 and dg(6,y)/ 
dy a /3 = &apl '(ra aa )Y,, where 8 a p is the Kronecker delta. Let So = (<T a 6,o) 
and Sq 1 = (oq*) be the true matrices of £ and X -1 , respectively. Then, 
since 5(A) = E /(A), dg{X)/dy a p = 5 af s/{ra aafi )E and g~ x (X) =E " 1 //(A), 
from (5), 

2/(A) 2 \r-CT Q! a,0O"77,0 ro- aQi0 rcr 77 ,o 
4 

= Ma^, j( A), say. 

i=l 



Thus, 



a/3 7 ^t 

where r = Ea,b=i ^o/K^o^.o) 
Next, from (6), 



\ fir j 

J2 ~ ^af3 1 y,i{X)g ap (X)g ll3 {X)dX = — , 
tt Jo Zr 



~ / / i a/37^2(A)ffa 1 /(A)g 7/3 (A)(iA 

a/37 At ° 
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~J Ma/37^3(A)5ai/(A)<7 7j 3 (A) d\ 



Finally, from (7), 



Similarly, 



1 

~27 



a/3 
1 

2' 



Consistent estimators of 9, r\ and a 2 are obtained by substituting S n 
(<5"a6,n) = h T,t=l Z * Z t for S - 



Example 2 (T/ie independence between the component time series). The 
independence between the component time series of a stationary Gaussian 
multivariate time series is equivalent to 

/ a6 (A) = 0, a / b, 

for all A and a, b = 1, 2, . . . , r. Thus, v = and g(y) = diag(yn, . . . , y rr )- 

Hence, g~ l {y) = diag(yf 1 1 , . . . ,y~ l ) and dg(y)/dy a p = 5 a p diag(5i a , . . . , 
S ra ). Then, since g(X) = diag(/n(A), . . . J rr (X)), dg(X)/dy af3 = 5 a p diag(<5i Q , 
...,S ra ) and g _1 (A) = diag(/f 1 1 (A), . . . , /~ 1 (A)), from (5), 

- 2 V/ Q q(A)/ / 3 /3 (A) " / QQ (A)2 J ■ 

Thus, 77 and cj 2 are the known constants expressed as 

V=j{r 2 -r), a 2 = l(r 2 -r). 

Hong [16] and Eichler [11] consider the same problem for the bivariate time 
series including non-Gaussian time series. Our result is a generalization to 
a multivariate one, though being derived under Gaussian assumption. 

Finally, we consider an example where a nonlinear constraint is imposed 
on the null hypothesis. 



12 



Y. YAJIMA AND Y. MATSUDA 



Example 3 (The conditional independence). Set Z a = {Z a t, —oo<t< 
oo},Y ab = {F a6jt ,-oo < t < 00} where Y abjt = {Z jt ,j / a,b} is an (r - 2)- 
dimensional random vector. Then, the conditional independence between 
Z a and Z b given Y ab is defined by 

(8) Cov(e a \ {atb}c {s),e b \ {a)b} c{t)) = 

for all s,t£Z = {0,±l,±2,...}, where 

00 

e a\{a,b} c (*) = Z at~ Y d* a (t - u)'Y abjU , 
u=—oo 

and d*(t — u) is the (r — 2) -dimensional vector which minimizes 

/ 00 \ 2 

E(z at - Y, d a {t-u)'Y ab A . 

\ u=—oo / 

Then, the relation (8) is equivalent to 

/ ab (A) = 0, -vr<A<^ 

(see, e.g., Dahlhaus [7]). 

Now, let V = {1, . . . , r} and £ be a subset of V x V. Consider the null 
hypothesis 

/ afe (A) = 0, (a,b)<£E 
for all A 6 (— vr,7r]. Then, u = and 

9ab{y)=yab, (a,b)eE, 

(9) 

<f fc (y)=0, (a, 6)^. 

A graph G = (V, £7) is the partial correlation graph of a time series Zt if 
(a,b) (£Eof ab (\)=0,y\ £ (-7r,vr] (see Dahlhaus [7]). It is proved by the 
implicit function theorem, in the same way as Lemma 7 of Matsuda, Yajima 
and Tong [26], that the constraint (9) is expressed in the form (1). 

Next, we show the outline of the derivations of n and u 2 , because it is 
given in the same way as Theorem 3 of Matsuda, Yajima and Tong [26] and 
Theorem 2 of Matsuda [24] by applying Lemma 8 of Matsuda, Yajima and 
Tong [26]. 

Let M = #{(a, b)\(a, b) ^ E,a <b} and 1q be the indicator function which 
is 1 if the condition C is true and otherwise. Then, we have 

(A) 



1 

2 7 (a,/3)e£^( 7 ^)e£; 

ab 



ka^eEh-r^E Y ( zZ ka,e)iE -JT <f & (A) + Ia=a/ b (A) 

. u yaj3 J ae , 



zZh 



l (b,e')<£E 
e> 



dgW 
dg(X) 



dy, 



g ea {\) + I b=1 g» a {\) 

be' j 
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Thus, 




=- M + -L r[tv(i r )] 2 d\ 

27r Jo 
= —M + -r 2 . 



From (6) 



r 



7T 



E 



A t a/3 7 ^2(A)^(A)5 7/ 3(A)(iA = - -(r 2 - 2M). 



a, /3, -y,u=l 



Then, it follows from (7) that 77 = M/2. Similarly, a 2 = M/3. 

Remark 2. All of the y) in Examples 1-3 map the space of positive 
definite matrices into itself (see Dempster [9] for Example 3). However, it 
is not necessary for g(0, y) to satisfy this condition under alternatives, be- 
cause Theorems 1 and 3 depend only on the behavior of g(9,f(\)) in the 
neighborhood of the null hypothesis, and Theorem 2 still holds by assigning a 
sufficiently large value to T n if fnj is not positive definite, and, consequently, 
K(M t ) cannot be defined. For practical use, being not positive definite, fuj 
gives strong evidence against the null hypothesis, and, consequently, we can 
reject it without causing any serious problem. 

5. Simulation results. We conduct some computational simulations to 
see the performance of the test statistics. Consider the following three- 
dimensional model for testing the independence between the component 
scries: 



where et are independent normal variables with mean and covariance ma- 
trix I3. The component series of (10) are mutually independent if (f> = 0.0 
and are dependent otherwise. 

Applying our test statistic T n with Kj(A) and u(x) = 1 on [—1/2, 1/2], we 
test the independence between the component series mentioned in Example 2 
of Section 4. We also considered Wahba's test statistic T* and the quadratic 
test statistic Tq^ for comparison. We examine performances of the tests 
under the null and alternative hypotheses in Tables 1 and 2, respectively. 



(10) 




0.7 (f> 0.0 
0.0 -0.5 <j) 
0.0 0.0 0.6 
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Table 1 shows mean, variance, skewness, kurtosis, 5% upper quantile of the 
null distributions and the empirical sizes under 5% asymptotic significance 
level based on the 1000 replications of the process (10) with = 0. We 
find that all the quantities except for skewness converge to the asymptotic 
limits relatively fast. However, it should be pointed out that bias exists 
in every test statistic in the small sample size, which causes nonnegligible 
size distortions. Hence, for fair power comparison, all the sizes of the test 
statistics are adjusted to be exactly 0.05 in the same way as Haug [15] 
and Saikkonen and Luukkonen [34]. Table 2 shows the empirical powers of 
the three statistics based on the 1000 replications of the process (10) with 
<fi = 0.1 and 0.2. We show the results when the bandwidth is predetermined 
as n = 101, m = 16 and n = 201, m = 30 and is selected by cross validation, 



Table 1 

Comparison of the null distribution of the tests for independence where 5% size is the 
empirical frequency of the rejection by the asymptotic 5% critical point 



Test 


Sample size 


rn 


Mean 


Var 


Skew 


Kurt 


5% qtl 


5% size 


f n 


101 


16 


0.163 


1.04 


0.66 


4.01 


1.94 


0.075 




201 


30 


0.097 


1.02 


0.62 


3.61 


1.94 


0.083 




501 


60 


0.073 


0.94 


0.59 


3.62 


1.74 


0.059 




1001 


120 


0.069 


0.94 


0.44 


3.09 


1.75 


0.064 




101 


16 


-0.163 


0.78 


0.59 


3.70 


1.44 


0.036 




201 


30 


-0.084 


0.87 


0.59 


3.56 


1.63 


0.048 




501 


60 


-0.028 


0.87 


0.59 


3.63 


1.60 


0.045 




1001 


120 


0.020 


0.90 


0.42 


3.06 


1.67 


0.053 


ft* 
± tt 


101 


16 


0.161 


1.14 


0.76 


3.91 


2.07 


0.095 




201 


30 


0.073 


1.00 


0.69 


3.43 


1.91 


0.075 




501 


60 


0.067 


1.01 


0.64 


3.41 


1.80 


0.070 




1001 


120 


0.050 


1.01 


0.47 


3.08 


1.89 


0.072 




oo 




0.0 


1.0 


0.0 


3.0 


1.64 


0.050 



Table 2 

Power comparisons of the size adjusted tests for independence under 5% 
empirical significance level 



Sample size 


4> 




Tq, n 


■*- n 














m = 16 






CVLL 




101 


0.1 


0.130 


0.118 


0.097 


0.110 


0.118 


0.098 




0.2 


0.503 


0.482 


0.275 


0.451 


0.451 


0.276 








m = 30 






CVLL 




201 


0.1 


0.231 


0.231 


0.163 


0.206 


0.232 


0.145 




0.2 


0.858 


0.859 


0.719 


0.831 


0.840 


0.621 
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which minimizes 

i [n/2] 

CVLL(m) = - ]T (WzJu)-i) +logdet(/ t/jW )), 
n 3=1 

where 

j m/2 

U,j,—j = — Yj J Z,j+k- 
k=-m/2,k^0 

Alternatives to CVLL can be the methods proposed by Lee [21] and Ombao 
et al. [28]. For tentative comparison between the performance of the fixed 
bandwidth selection and the data driven one, we adopt CVLL. However, we 
should remark that CVLL can select the bandwidth minimizing the mean 
squared error asymptotically if we put f3 = 4/5, which lies outside the interval 
of (A4) (see, e.g., Beltrao and Bloomfield [1], Matsuda and Yajima [25]). This 
issue often emerges in nonparametric hypothesis testing (see, e.g., Zhang [41] 
and the references therein and Fan and Yao [12], Section 9.2.7) and more 
rigorous consideration is left to future studies. 

We find from Table 2 that our statistic has almost the same power as 
the quadratic one, whereas the Wahba's is significantly less powerful than 
the others, which reinforces that Wahba's statistic is less efficient than 
the others in Pitman's ARE sense as mentioned in Remark l(iii), because 
2D U /C U = 2/3 is less than 1. 



6. Proofs of theorems. 



Proof of Theorem 1. Applying the Taylor expansion and Lemma 2 
to K{Mt) and noting that the first-order terms of the expansion vanish, we 
obtain 

1 [n/2] 
L i=l 



(11) 



+ g 2^ K ab,cd,ef 
a,b,c,d,e,f=l 

[n/2] 



x (m a b,t ~ <5afe)("i c d,i - Scd)(m e f,t ~ <M 



+ Op(l). 



It is shown later that the second term on the right-hand side of (11) is o p (l). 
Hence, we shall consider the limiting distribution of the first term. 



16 Y. YAJIMA AND Y. MATSUDA 

Note that /(A) = 5(A) if the null hypothesis (1) is true. Then, it follows 
from Lemma 2 and 9 n — 9q = O p {n~ l l 2 ), by the Taylor expansion, that 

M t - I r 

= fuAfnb ~ 9t l ) + (fu,t - gtjgT 1 
= gt{g- 1 (O n Ju,t)-g- l (0 j t )} 
+ (fu,t - gt){g'\o n J U)t ) - g -\e Q j t )} + (f Ujt - gt )g^ 

12 

r v 

= 3-tidg.t 1 /dy a p){f U)af 3t- fa/3,t) + ^2gt(dgr 1 /d9i)(6 in - e i0 ) 

a,/3=l i=l 
1 r 

+ 2 Y 9t{d 2 gt l /dy a pdy lv )(f Uta pj - f a /3,t)(fu,<yv,t - fju,t) 

+ 0p ( n -3/4) + { f ut _ g t ){g-\e n J U)t ) - 9 ~ l (9 , f t )} 

+ (fu,t- gt)9t X - 

By applying Lemmas 2, 6 and 7 and 9 n — 9q = Opin- 1 / 2 ) to (12), 
£tr[(M t -/ r ) 2 ] 



(13) 



Ydt( d 9t 1 /dy af 3)(fu,a(3,t - faP,t) + (fu,t ~ 9t)9t 1 ) 

/ . 



+ o p (n 1 / 2 /m 1 / 2 ) 

say. 

Next, define by replacing f ut - g t and f Utt - f t of Tl 1] with ft - E{ft) of 
Section 7. Then, from Lemma 8, it suffices to derive the limiting distribution 

of \TP. 

By an elementary calculation, 

1 r 

where H a p-yu,t,n = Va(3yv(^t)- Then, if we define f n by 

V n \2 m J 

it follows from Lemmas 9, 10 and 11 that lim n E(T n ) = 0, lim n Var(T n ) = a 2 
and for k > 3, the kth order cumulant of T n converges to as n — > oo. Hence, 
the limiting distribution of T n is iV(0,cr 2 ). 
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Finally, by the same argument as (12), it is shown that the second term 
on the right-hand side of (11) is negligible. □ 

Proof of Theorem 2. f n is written as 



'nm 



Vn <r n - 

n m J 



Then, it follows from Lemma 2 and Exercise 1.7.4 of Brillinger [2] that 
T n j (na n ) converges in probability to Jq K(f(X)g(9* , /(A)) -1 ) dX as n — > 
oo. Thus, the assertion is shown immediately. □ 

PROOF of Theorem 3. First, we note that, under H an , Lemmas 1-11 
of Section 7 are still true if we substitute ft, n for ft in the definitions of yt 
and y\ and add the term of 0(l/(mn) 1 / 4 ) to the right-hand side of Lemma 
9. Next, since M t — I r = o p (?i -1 / 4 ) is also true under H an , analogous to (13), 
we have 

]>>[(M t - I r f] 



+E tr 



^ 9t,n (p9t.n I dVap) (fu,a0,t ~ fa(3,t,n) + (fu,t ~ ft,n)9t,n 
L\a(3 / 

2-1 



(ft - gt,n)g+n + 



1 



(nm) 1 / 



\u ft 9t,n 



+ o p (n 1 / 2 /m l l 2 ) 



say. 



Then, by the Taylor expansion and Exercise 1.7.14 of Brillinger [2], 



n 



1/2 



Hence, 



(14) 



m 



n 



1/2 



(20+o(n 1 / 2 /™ 1/2 )- 



T n ~ — fin) I O n 

m ) 
m /l-(i) n \ 
n \2 ml 



+ £/&n + \ — {(V ~ Vn) + (Vn ~ Vn)}/^n + O p (l). 

V m 

The limiting distribution of the first term on the right-hand side of (14) is 
the standard normal distribution. Then, the assertion is obtained, because 
r/n -r ? = 0(l/(mn) 1 / 4 ) and o* - a 2 = 0(l/(mn) 1 / 4 ). □ 
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7. Lemmas. First, we introduce some random variables, notation and 
technical remarks. Define y t = {y a b,t),Vt = (y a b,t),t = l,...,[n/2], by y t = 

fu,t ~ ft and y t = fu,t - E(fu,t), respectively. 

(A2) and (A3) assure that J^ 7r logdet/(A) dX > — oo. Then, under (Al), 
Z t is expressed as 

oo 
3=0 

where $j is an r x r matrix and J^jLo^i^j^'j) < °° and £t = (eit, • ■ • , £ r t)' 
is a mutually independent zero-mean Gaussian process with covariance ma- 
trix I r . Then, denote the discrete Fourier transform of e a t by W £ {\) = 
^"=i £at x ex P(^A) for a= l,...,r and the cross periodogram of e at 

and eu and the periodogram matrix by ^(A) = W £ (X)W b (X) and I e (A) = 
{ ie abW)i respectively. 

Next, define the r x r matrix ff = {f £ bt ) by 

1 m/2 

fab,t = — ^^a(exp(a t+j ))/f +j ^(exp(-iA t+j )), l<t<[n/2], 

j=-m/2 

where Jj = / e (Aj) and $ a (e a ) is the ath row vector of <£>(e iA ) = E£L %e^ A . 

Then, define ?/f = / t £ — ft and = ft — E(ff), respectively. We shall 
derive the main results by showing that the limiting behaviors of T n remain 
unchanged if we substitute yfiVt) f° r Vtivt)- Vt aim Vt are more tractable 
than y t and yt, because 2irlj ,j = 0, 1, . . . , [n/2], are mutually independent 
random variables with Wishart distributions and, consequently, y\ (yf) and 
y £ s {y £ s ) are mutually independent for \t — s\ > m. 

Lemma 1. Under (A3) ; 

E(W ad W bJ ) - f absj = O , -m/2 + 1 < j < [n/2] + m/2, 

E(W a>j W b:k ) = O i -2- J , -m/2 + 1 < j < k < [n/2] + m/2, 

j + k^0,n, 

/log n\ 

E(W aJ W b ,k) ~ fab,j = O MM - -m/2 + 1 < j < k < [n/2] + m/2, 

j + k = 0,n, 

E(WajW blk ) = O i -M , -m/2 + 1 < j < fc < [n/2] + m/2, 
uniformly in j and k for a, b = 1, 2, . . . , r. 
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Proof. The assertion is shown in the same way as Lemma 1 of Matsuda 
and Yajima [25] by noting that W a ,jW bt k = W a ,jW bi j for j + k = 0,n. □ 

Lemma 2. Under (Al)-(A4), 

sup \fu,ab,t ~ fab,t\ = o p (n" 1/4 ). 

i=±l,...,±[n/2] 

Proof. Applying Lemma 1, we see that Lemmas 2 and 3 of Matsuda 
and Yajima [25] are still true for fu,ab.t arid fabf Then, the assertion is 
proved in the same way as Proposition 1 of Matsuda and Yajima [25]. □ 

Lemma 3. Under (A3), 

E(y ab>t ) = 0(m 2 /n 2 ), 
uniformly in t = 1,2, ... , [n/2], for a, b = 1, 2, . . . , r. 

Proof. It follows from Lemma 1, by the Taylor expansion, that 
1 m/2 

E(y a b,t) = — u 'j(fab,t+j + O (log n/n)) - f a b,t 

W i=-m/2 

1 ^ / 27T7' \ 

— E W A f»b,t + —fkt + 0(f/n 2 ) + 0(logn/n) - f ab , t 



3 =-m/2 

= 0(m 2 /n 2 ) + 0(logn/n) 
= 0(m 2 /n 2 ), 
uniformly in t where f' ah t = df ab (\)/d\\ x =\ t . □ 

Lemma 4. Under (Al), (A3) and (A4), 

, , ( 0(l/m), if\t — s\<m, 

Cov(y aM , y cc( , s ) - | 0(lQg2 n/n2)) . f |t _ s| > ^ 

uniformly in t, s = 1,2, ... , [n/2], /or a, 6, c, d = 1, 2, . . . , r. 

Proof. We observe that 
Cov(y ab ,t,ycd,s) 

= E(y abtt y cd:S ) - E(y abit )E(y cd)S ) 

1 



m/2 



E ^Wfc[^(^a,i+j^ C , S +fc)£'(W / M+jW / " (iiS+ fc) 
V" 7 j,k=-m/2 
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Then, the assertion follows immediately from Lemma 1 by noting that t+j 7^ 
s + k and t + j + s + k 7^ 0, n for any j and k if |t — s| > m. □ 

Lemma 5. f/nder (Al), (A3) and (A4), 

cum(y ai6litl , y a2 62,t2 > • • • > Va k b k ,t k ) = 0{m l ~ k ), 
uniformly in 1 < ti < [n/2] for any k>3 and 1 < otj, 6j < r, z = 1, . . . , fc. 

Proof. The assertion is obtained, in the same way as the proof of The- 
orem 7.4.4 of Brillinger [2], by applying Lemma 1 and Theorems 2.3.1 and 
2.3.2 of [2]. □ 

From now on, let /x t)Tl , t = 1, . . . , [n/2],n = 1,2, ... , be constants bounded 
in t and n. 



Lemma 6. Under (Al), (A3) and (A4), 

[n/2] 

fh, n yab,t = O p (n 1/2 ) 

t=l 

for a, b = 1, . . . , r. 



Proof. We may assume, without loss of generality, that ^t,n = 1. Then, 
it follows from Lemmas 3 and 4 that 



E 



J2Cov(y ab:t ,y ab:S ) + 



s.t 



E(Vab,t 



= 0(n) + 0(log 2 n) + 0{m A /n 2 ) 
= 0(n), 

which implies the assertion. □ 

Lemma 7. Under (Al), (A3) and (A4), 

[n/2] 

Y fH,nl/ab,tycd,tyef,t = Op(l) 
t=l 



for a,b,c, d, e, f = 1, . . . , r. 



Proof. We may assume, without loss of generality, that fJit,n = 

1. It 

follows from Lemmas 2 and 3 that 

y a b,t = Vab,t + 0(m 2 /n 2 ) = o p {n~ l/A ), 
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uniformly in t. Thus, 

(15) Vab,tVcd,tVeUt = X] Vab,tVcd0ef,t + O p (l). 

t t 

Now, we evaluate the first term on the right-hand side of (15). We have 

2 

Vab,tVcd,tVef,t 



E 



(16) 



X Cov(y abtt ycd,tyef,t, yab, S ycd, S yef,s) + 



X E(y a b,tycd,ty e f,t) 



Noting E(y abt ) = and applying Theorems 2.3.1 and 2.3.2 of Brillinger [2], 
we observe that the first term on the right-hand side of (16) is equal to 

X cum(y ab jy cdtt y ef j, y b a, s ydc,syfe,s) 
s,t 

= ^[cum(y afeji , y cdjt , y ef>t , yba,s,ydc,s,yfe,s) 
s.t 



(17) 



+ cum(y afe> i , y cd j) cum(y e / >t , y ba , s ,ydc,s, Vf e ,s) 
+ cum(y aM , y ba:S ) cum(y cd:t ,y ef j, y dCtS ,yf e ,s) 

+ cum(y ab j,y c d,t, yba,s) cum(y e / )t , y dc , s , yf e>s ) 

+ cum(y aM , y cd:t ) cum(y 6ajS , y dCtS ) cum(y eftt ,y f ^ s ) 

+ cum(ju i4 , y ba:S ) cum(y cdjf , y dCtS ) cum(y e/]t , y /6)S ) 

+ the remainder terms] . 

Any term of the remainder ones on the right-hand side of (17) is expressed 
in the same form as one of the preceding six terms. Hence, it suffices to 
evaluate these terms. It follows from Lemmas 4 and 5 that the first one 
is 0(n 2 /m 5 ), the second and fourth ones are 0(n 2 /m 4 ), the third one is 
0(n/m 3 ) and the fifth and sixth ones are 0(n/m 2 ), respectively. 
Similarly, the second term is equal to 



X] cum (y~ab,t, ycd,t, y e f,i 



^cum(y ab j,y cd j,y eU ) 



= 0(n 2 /m 4 ). 

Thus, the proof is completed. □ 
Lemma 8. Under (Al)-(A4), 
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Proof. It suffices to show that 

[n/2] 

(!8) Y teAVaPjahn*,* ~ y E a/3,ty*v,t) = o p {n l / 2 /m 1 / 2 ) 

t=i 

for a,j3,-),v = l,...,r, because — T n 2 ' is a summation of the terms 
expressed in the same form as (18). We may assume that Ht n = 1 without 
loss of generality. By the Schwarz inequality, the absolute value of the term 
on the left-hand side of (18) is bounded by 

1/2 / \ 1/2 



(Y \VaP,t ~ fal3,t\^j (Y \y^t\ 2 ^j 

+ (Y\y £ afs, t \^j (zk- 



(19) 

\ 1/2 

The first term of (19) is bounded by 

(20) (^Ylya^t-y^f + ^Yly^t-fa^ (yM 2 ^ ■ 

It follows, from Lemmas 3 and 4 and the assertion during the proof of 
Lemma 2 of Matsuda and Yajima [25], that the first term of (20) is equal 
to O p ( (log n/m) 1 / 2 ) + 0(m 2 /n 3 / 2 ) and the second term is O p ((n/m) 1 / 2 ). 
The second term of (19) is evaluated in the same way. Thus, the proof is 
complete. □ 

Lemma 9. Let n(X) be a differentiable function in [0,7r] , and set nt, n = 
fi(X t ),t=l,...,[n/2]. 

Then, under (Al)-(A4), 

m [n/2] c u r 

— Y E (Pt,nfab,ty £ cd,t) = — M(A)/od(A)/ c6 (A) dX + 0(m/n) 
n ' it Jo 

for a, b, c, d = 1, . . . , r where 

1 fl/2 / fl/2 \ 2 

C u = - u(x) dx/[ / u(x)dx) . 

2 J -1/2 \J-l/2 ) 

Proof. Applying Lemma 1 to W^j,a = l,...,r, and noting that the 
right-hand side terms of the lemma are exactly equal to for W^j, we 
obtain 

[n/2] 

— Y E ^ny £ a h , t y £ cd,t) 
u 1=1 
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m x 2 [n/2] m/2 

ZJ^v (9^) ^ ^ 2 ^,n$a(exp(At +i ))^(exp(-A t+J )) 
j \ z7T / t=1 j= _ m/2 

x $ c (exp(A t+i ))$ / fe (exp(-Ai +i )) + 0(m/n) 



[n/2] m/2 

H w "jVt,nfad,t+jfcb,t+j + 0{m/n) 



n(w*) 2 ~^ 



m 



j=-m/2 
m/2 [n/2] 



Y w ] Y Vt,nfad,tfcb,t + 0(m/n). 



n(w*) 2 . 

The last equality is given by the Taylor expansion. Then, the assertion fol- 
lows from Exercise 1.7.14 of Brillinger [2]. □ 

Lemma 10. Let /ii(A),i = 1,2, be differentiable functions in [0,n] and 
set fj,i t t, n = /^i(At) . Then, under (A1)-(A4), 

m /[n/2] [n/2] x 

\t=l 8=1 / 

D ± r W (A)MA) 



vr Jo 

X (/oia 2 (tyfb 2 bi (A)/cic 2 (tyfd 2 di (A) 

+ /aic 2 Wfd 2 bj (A)/ Cl a 2 (A)/&2di (-^)) ^ 
/or ai,bi,Ci,di = 1, . . . , r, i = 1, 2, uTiere 

! /-1/2 ,1/2 ,1 / ,1/2 \4 

Du = t; \ dx dy I u{x)u(y)u{x + z)u{y + z) dz / \ I u(x)dx) . 
2 7-1/2 7-i/2 7-i V7-1/2 / 

Proof. First, note that 4 and y e cd s are mutually independent for \t — 
s\ > m, because ij, j = 0, 1, . . . , [n/2], are mutually independent variables. 
Then, noting that E(y £ ab t ) = and applying Lemmas 1 and 5 and Theorems 
2.3.1 and 2.3.2 of Brillinger [2], we obtain 

^ Cov (Y Vl^nVliwMidi ,t> Y^^fa^sfc^^ 

J2J2 CUm ^t,nfa 1 b 1 Mld u uV^M 2 a2,8yd2C2,s) 



m 

71 t s 

m 
n 



t s : \t— s\<m 
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+ cum(y £ aibut , f b2a ^ s ) cum{tf Cldut , #| 2C2 , S ) 
+ cum(^ l6ljt , y £ d2C2)S ) cum(y £ Cidut f b2a2 S ) 

Tfl ^ — > ^ — > ^ — -v 



n(w*) 4 ~ . , . r-r* . 

v ' t S,\t-s\<mn J2J3J4 

X i.faia2,tfb2bi,tfciC2,tfd2di,t 
X ^+il,S+j 3 ^+j2,S+i4 
H" faiC2,tfd2bi,tfcia,2,tfb2di,t 

x ^+ii,s+i4^+i2,s+is) + (!) 

* |t— s|<m jij2 

X (.faia,2,tfb2bi,tfciC2,tfd2di,t 

+ fa 1 C2,tfd 2 b 1 ,tfc 1 a2,tfb 2 d 1 ,t) + o(l). 

The last two equalities follow from 2tt ji/n, i = 1, . . . ,4 and 2ir(t — s) /n being 
of order m/ra, and, in the last summation, ji,j2,s and t must satisfy +£ — 
s| < m/2 and |j2 + i — s\ < m/2. Then, the assertion is obtained by Exercise 
1.7.14 of Brillinger [2]. □ 

Lemma 11. Set 

T abcd,n — I I 2-~i c obcd,t ) nyab,tycd,t 

for a, b, c, d = 1, . . . , r, where c a b c d,t,n are the constants bounded in t and n. 
Then, under (Al)-(A4), 

cum ( T a 11 bna 1 2bi2,m ■ ■ ■ i T a kl b kl a k 2b k2 ,n) = 

for k>3 and 1 < a,ij,bij < r, i = 1, . . . , k, j = 1, 2. 

Proof. Applying Theorems 2.3.1 and 2.3.2 of Brillinger [2], we have 

CUm(r ail ;, liai2 b 12jn , . . . , T a kl b kl a k2 b k2 ,n) 

) 11 c a il b il ai2b i 2,t i ,n 

J h,t 2 ,...,t k =li=l 

( 21 ) X CUm (^a 11 b 11 ,t i ya 1 26i2,ii'---'ya fel fe fel ,t fe 2/a fe2 6 fc2 ,t fe ' 



\ fc/2 H 2 1 k 

e n 

Caabiiai2bi2,U,n 

7 h,t 2 ,...,t k =li=l 
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v 
Z=l 

where 2~^* is over au the indecomposable partitions v* = v\ U • ■ • U u p of the 
table 

11 12 

21 22 

fcl jfe2. 

Let I* be the number of elements of vi, I = 1, . . . ,p. Since E(y £ abt ) = 0, any 
/* is greater than 1, and, hence, p satisfies p< k. 

Now, applying Lemmas 4 and 5 to each cumulant on the right-hand side 
of (21) and noting that J2f=i = 2fc, we have 

f[cum(f ai . bijtti ,ij G ^) =o(nm w *) = 0(m^ 2fc ). 
z=i \z=i / 

By Theorem 2.3.1 (iii) of Brillinger [2], the number of the nonzero terms of 
the summation T,t u t 2 ,...,t k is 0(nm fc " 1 ), since //„ ,, ., and . lh ., are 
mutually independent for \ti — ti>\ > m. Thus, the term on the right-hand 
side of (21) is 0(mP~ k / 2 ~ 1 /n k / 2 ^), which is 0((m/n) k ^ 1 ) = o(l) for k > 3. 
Hence, the proof is complete. □ 
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